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Abstract 

The cross section through gluon fusion is calculated for the production of the light 
neutral Higgs boson through next-to-leading order QCD within the Minimal Super- 
symmetric Standard Model. The quark-mediated contributions are taken into account 
exactly, while for the genuinely supersymmetric terms we use expressions obtained in 
the limit of large squark, gluino and top quark masses. We present numerical results 
for the total inclusive cross section as well as for kinematical distributions of the Higgs 
boson. We also consider the effect of an MSSM-like 4**^ generation on the total Higgs 
production cross section. 

1 Introduction 

With the first sensitivity of the Tevatron to a Higgs signal and the first data taken 
at the Large Hadron Collider (LHC), Higgs phenomenology has received an additional 
significant boost in the last few years. It is pleasing to see that the efforts towards precise 
predictions for Higgs production (for reviews, see Refs. [IHS]) at the LHC are paying off, 
currently leading to significant bounds on the Standard Model (SM) Higgs boson mass. 

Another interesting application of the on-going Higgs searches are the restrictions on pa- 
rameters of theories that go beyond the SM, such as the Minimal Supersymmetric Standard 
Model (MSSM) [6] (see also Ref. [7]) or the SM with a fourth fermion generation [8].^ 

One of the most important production mechanisms for these searches is gluon fusion 
which is known in the SM through next-to-next-to-leading order (NNLO) QCD [9^119]. 
Effects beyond that order have been investigated as well [20H25] . with the reassuring 
conclusion that the theory uncertainty derived from the renormalization and factorization 
scale variation of the fixed order NNLO result seems to be reliable. Electro-weak (EW) 
corrections have been found to be below around 6% relative to the leading order (LO) 



^ See |http : //tevnphwg . f nal . govT] for updates. 
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result |261I27|. and a calculation of the mixed QCD/EW effects in the unphysical limit 
Mfi <C Mw supports the assumption of a factorization of QCD and EW corrections [28j . 

Concerning the MSSM, at LO squark loop effects have to be taken into account. Due to 
supersymmetry (SUSY), their couplings to the Higgs bosons are typically suppressed by 
powers of Mg/Mq, which is why they only contribute significantly to the cross section if 
the squark mass Mq is not too large. In that case, however, there are certain regions of 
the SUSY parameter space where a quite drastic cancellation among the quark and the 
squark induced amplitudes may occur |29ll30j . 

Apart from the SM-like diagrams involving only quarks and gluons, QCD corrections to 
gluon fusion in the MSSM require the calculation of the corresponding corrections to the 
LO squark diagrams, plus diagrams involving quarks, squarks, and gluinos simultaneously. 
Meanwhile, all the ingredients for the next-to-leading order (NLO) QCD corrections to 
gluon fusion in the MSSM have been calculated: The top/stop/gluino effects were taken 
into account through an effective Lagrangian approach |31H34] . a result for scalar QCD 
(scalar quarks with arbitrary mass) was calculated both numerically and analytically [351 - 
I37j . and for the bottom/sbottom/gluino diagrams recently a result in terms of an expansion 
in the inverse SUSY masses was presented [38]. The quark/squark/gluino effects have also 
been calculated fully numerically both for the top and the bottom sector [39]. Very 
recently, even the NNLO effects in the heavy top/stop limit have been evaluated |40j . 

In this paper, we present a calculation of the full NLO QCD effects for hadronic Higgs 
production in the MSSM with real parameters and combine it with the NNLO QCD effects 
from purely top-quark induced contributions known from the SM calculations. Although 
most of our results should be valid for Higgs masses up to < 2min(Mf, M) (M is the 
typical scale of the SUSY partner masses), we present results only for the phenomenolog- 
ically most relevant case of the light MSSM Higgs boson, for which Mh < 135 GeV (see, 
e.g. Refs. [41^ 142) for reviews, and Ref. [l3] for the most up-to-date result). We have 
re-calculated all of the real corrections due to quarks and squarks, while for the virtual 
corrections to the pure quark terms the analytic result of Ref. [H] is employed. These 
results are valid for general quark/squark masses and can therefore be used for the top 
as well as the bottom sector. Concerning the virtual squark effects, we use the effective 
Lagrangian approach of Ref. [32] for the top sector, evaluating the Wilson coefficient with 
the help of evalcsusy . f [35] . For the bottom sector, we follow the approach of Ref. |38II46| 
in order to evaluate pure sbottom and the bottom/sbottom/gluino diagrams in terms for 
l/Mfo, where Mf, denotes the average sbottom and gluino mass. 

Apart from the inclusive total cross section, we also present kinematical distributions for 
the Higgs boson: while the transverse momentum (px) and pseudo-rapidity [r]) distribu- 
tions at 0{ag) have to be considered as LO (at 0{a'^), it is pT = and r] = oo) and 
have been evaluated before [UHlH], to our knowledge this is the first time that rapidity 
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(y) distributions are being presented in the full MSSM at this order. Our calculations are 
based on a numerical routine which will be made publicly available and can be obtained 
from the authors upon request. 

The remainder of this paper is organized as follows: Section [2] outlines our notation, 
describes our calculation and methods, and quotes some of the most important formulas. 
Section [3] defines the scenarios for our numerical analyses, describes how we obtain our 
best prediction of the cross section, and provides numerical results for the total inclusive 
cross section. Section H] presents results for the LO transverse momentum and the NLO 
rapidity distribution of the Higgs boson. Section [5] considers the effect of an MSSM- like 4**^ 
generation on the total inclusive cross section through NLO QCD, and Section [U] contains 
our conclusions. In the Appendix, we collect the Feynman rules for the Higgs couplings 
used in this paper, give some more analytical formulas, and provide numerical results for 
other SUSY scenarios. 



2 Calculation of the cross section 
2.1 General outline 

For the sake of clarity, we collect some of the most important notation at this point. 
Concerning the SUSY parameters, we follow the usual notation (as mentioned before, we 
work in the MSSM with real parameters): 

• The angle a. rotates the CP even neutral components of the two Higgs doublets }i\ 
and into their mass eigenstates /i, -ff, where by definition < Mh- In all our 
formulas, a actually denotes the loop-corrected effective mixing angle a^s- 

• The ratio of the two Higgs vacuum expectation values is denoted by tan/3 = V2/V1. 

• The coefficient of the bilinear term ~ Hi ■ H2 in the superpotential is called ^SUSY- 
At tree- level, it is related to other SUSY parameters as in Eq. (p6]l of Appendix Rl 

Furthermore, the following abbreviations will be applied: 

(1) 



AMI ~ -fS 
Ml ' ^^'^ Ml 



where Mq and Mqi denote the on-shell quark and squark mass (g S t}), while Mq = 
{Mqi + Mq2 + Mg)/3 is the average of the g-squark and gluino (Mg) masses. The mass 
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of the light CP-even Higgs boson is cahed M^, and fi^ is the renormahzation scale. The 
factorization scale will be denoted by /if. 

Unless stated otherwise, the strong coupling constant will always be defined in the frame- 
work of standard five-flavor QCD, renormalized in the MS scheme: 

= (2) 



Following the notation of Ref. jl2j . we write the hadronic cross section at the hadronic 
center-of-mass energy s through NLO QCD as 



a{pp H + X) = ao 



1 + 

vr 



Th-rz- + ^^gg + ^^oa + Act, 



'9Q 



qq ) 



where 



1 dx 



g{x)g{T/x) , 



(3) 



(4) 



with the gluon density g(x) and r/^ = M'^/s. The normalization factor ao determines the 
LO cross section. In the framework of the MSSM, we write it as 



288^/27r 



qG{t,b} 



(5) 



where 



4°^=5,^(l + (l-r,)/(T,)) 



;(0) 
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Y.9q,^^i^-rq,^f{fg,^)) , (6) 
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The a^^\Tq) and the d^^\fq^i) terms are due to quark and squark diagrams, respectively, 
as the ones displayed in Fig.[TJ It may be useful to quote the relevant limits of these 
functions: 
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(a) (b) (c) 

Figure 1: Feynman diagrams contributing to Higgs production in gluon fusion at LO. 





(a) 



(b) 



(c) 



Figure 2: Sample Feynman diagrams contributing to the real corrections at NLO to Higgs 
production in gluon fusion. The particle in the loop can be a quark or a squark of bottom 
or top flavor. 



The coupling constants and g^- are given in Appendix El those for q = t can also be 
found in Ref. [32]. 

In principle, the sum in Eq. ([5]) should run over all quark flavors. However, only the top 
and the bottom quark account for a relevant contribution to the cross section, while all 
others are suppressed by their Yukawa coupling. 

The quantities Acr^g, Acr^g, Aaqg in Eq. ([3]) denote the non-singular terms of the cross 
section arising from gg, gq and qq scattering. Typical diagrams are shown in Fig. [21 the qq 
contribution is obtained through crossing from Fig.[2](c). They can be evaluated using well- 
known techniques. We have expressed them in terms of Passarino-Veltman functions |50j 
and checked our result against the literature (see, e.g. Ref. [5THM] ). 



2.2 Virtual corrections 



The coefficient C in Eq. ([3]) denotes the virtual corrections to the gg initiated process, 
regularized by the infrared singular part of the cross section for real gluon emission. We 
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(a) 



(b) 



(c) 



Figure 3: Sample diagrams for the mixed quark/squark/gluino contribution to gluon fusion 
at NLO. 



write 

C = 2Re 



„2 

+ 7r2 + /3o In^^, (9) 
/if 



where /3o = 11/2 — n^/S with n/ = 5, and is the LO amphtude in the hmit of large 
stop and sbottom masses, i.e., 

The NLO quark-loop contribution corresponding to the SM part, has first been eval- 



uated numerically a long time ago [12] and was later expressed in terms of analytic func- 
tions [3511361H1]- We provide a few terms of its small- and large- mass expansions in 
Appendix [BJ 

One class of diagrams contributing to the SUSY part d'q''^ is obtained by attaching a virtual 
gluon to Fig.[T](b) and (c). However, both squark and gluino effects need to be considered 
in order to preserve supersymmetry |31ll32j : sample diagrams containing a gluino are 
shown in Fig.O A fully numerical result for a^^ for general quark/squark/gluino masses 
was obtained in Ref. [391, but the corresponding code is not publicly available. For the 
pure squark diagrams, there exists both an analytic and a numerical result |35H37| . In 
Ref. j3 H [32 l [34]. was evaluated for the top sector (i.e., q = t) in terms of an effective 
Lagrangian, and one can write 

dW=4^)-,,^^ + 0(M,^), (11) 

where c^^ is the NLO term of the Wilson coefficient, defined in Eq. (2.5) of Ref. [32]. It 
can be evaluated with the help of the publicly available program evalcsusy.f 
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Concerning the SUSY bottom quark/squark sector, a fully numerical result for general 
masses was presented in Ref. [39]. Recently, this contribution was calculated with the 
help of asymptotic expansions in the limit of large SUSY masses [38]. We present an 
independent result, following the same strategy. However, due to the SUSY mass spectrum 
emerging from the most popular SUSY breaking scenarios, we decided to take the limit 
Mb = Mbi = Mb2 = Mg which leads to an extremely simple result [Mg is the gluino mass). 
Let us briefly describe the calculation. 

For the pure sbottom diagrams (i.e., without gluinos), the procedure is very similar to 
the one applied in Ref. [T7HT9l [55 t [56] , where top quark mass suppressed terms for the SM 
cross section a{gg — )• H + X) were evaluated (mind you, through NNLO; here, we consider 
NLO only). The expansion of the mixed bottom/sbottom/gluino diagrams is a little more 
involved, but still rather straightforward due to the algorithmically defined method of 
asymptotic expansions (see, e.g. Ref. [57] )• Let us consider an example: 




The notation is as follows: in the original diagram, Taylor-expand all the propagators 
of the sub-diagram 7 left of ® in terms of p/ M before integration. Here, p denotes any 
dimensional quantity (mass or momentum) of 7 except its loop momentum or M. The 
resulting Feynman integrals to evaluate in the above case are therefore: (i) one- and two- 
loop tadpole integrals (i.e., vanishing external momenta) depending on M, but not on M^] 
(ii) one- loop vertex integrals with external momenta 9i = ~ and 2qi ■ q2 = M|, and 
mass Mf) (no dependence on M). Both types of integrals can be calculated analytically: 
type (i) with the help of MATAD [58], type (ii) by standard Passarino-Veltman reduction [50] . 
for example. 
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For clarity, let us consider another example: 



^ AAAA/ 




Using this procedure, the virtual corrections are valid for arbitrary values of Mh and Mf,, 
as long as they are both smaller tharg 2Mb. 

We perform the calculation in Dimensional Reduction (DRED) by explicitely taking into 
account e-scalars as propagating particles. In order to avoid infra-red singularities in the 
corresponding Feynman integrals, we assign a large mass to the e-scalars for which we 
can assume — )■ oo at the end of the calculation. Practical details about the implemen- 
tation of e-scalars, in particular their Feynman rules, have been given in Refs. |59fl61] . 

Concerning renormalization, we use on-shell conditions for the bottom quark mass, one 
of the sbottom masses, and the sbottom quark mixing angle. The other sbottom mass 
is then fixed by SU(2) symmetry. The gluino does not require renormalization at this 
order. This top/stop-like renormalization scheme may cause perturbative problems for 
large ^SUSY and tan/3 [62|. However, due to our limit of degenerate sbottom masses we 
will restrict our analysis to moderate values of tan /? anyway. The comparison of various 
other renormalization schemes is beyond the scope of the current paper and will be deferred 
to a forthcoming publication. 

The external gluon wave functions are renormalized on-shell using 




■^The factor of two can be deduced from the analytic structure of the amphtude. 
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where = 6 is the number of squark flavors, Ca = 3 and T = 1/2 are color factors, and 
the last term in 6z is due to the e-scalars. In fact, at this order of as, all of the finite part 
of Z3 cancels against the conversion (see, e.g. Ref. |63j ) 

= (1 - fc)afS,(5) 

where a^^ denotes the strong coupling constant renormalized in the full MSSM according 
to DR subtraction. 

The bottom mass Mh, the sbottom mass Mhi, and the sbottom mixing angle 9b are renor- 
malized on-shell in analogy to the top quark sector, see, e.g. Ref. [32]. The sbottom 
mass Mf,2 is determined from the SU(2) relation from the other (top and bottom) on-shell 
quark and squark masses and mixing angles, and its on-shell value M^^. The two mass 
definitions are related by a finite shift (see, e.g. Ref. [64j^ : 

M.2 = M.? (1 + f A^W) . (15) 



The sbottom sector contribution to the Higgs production cross section can now be written 
as {snb = sm{n0b), Cnb = cos{ndb)) 



sin a 

M2 cos /3 



MSUSY cos(a - /3) 
cos^ /3 



+ -4; sm a + /? 
Ml 



1 



Mb 

2 
— ( 
3 



OS 



■Sfo + 9^26 sin^ Gw ) 1^0^ 



(16) 



+ o(Mr^) 



with the weak mixing angle ^VK) cf. Eq. (j34p . Our result for the coefficients K\ is (cf. Eq. ([T])) 



25 17 



24 



jfin 



^'=24+12^'^^-- (^™^-^^^)+^ + "^(--^) 



17 r„ 



24 



'q) ) 



-^[l + (l-r,)/(T,)] + 



36f„ 



15. 



-37 + -l^s - 3r, - 3(2 - r, - r,^)/(r,) 



(17) 



Note that the /C? are given for a general quark fiavor q, so we will be able to apply them 
also for the top contribution further below. The only difference arises from the functions 



^Ref. fixes on-shell and evaluates M51 from the SU(2) condition; we found better numerical 
stability by fixing Mh\ instead. 
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/(r), defined in Eq. ([7]), and -Bo'^(r), given by 



rq 



2 — 2yjTq — 1 arctan 



- 1 



Tg > 1, 



«Vr , Tq < 1 . 



(18) 



Due to our limit of a degenerate SUSY mass spectrum, all the dependence on the squark 
mixing angle 6q drops out in Eq. (|17p . This also removes all terms ~ Mg/Mq found in Ref. 
|38| . In fact, many of the terms evaluated in that reference vanish in our approximation, 
while many of the terms of Eq. (|17p vanish in the approximation of Ref. [3S] since they 
are of higher orders in Mi, or l/M^. Needless to say that for the terms that are non-zero 
in both approximations we find complete agreement 

Since the bottom quark mass is very small, the following expansion in M^/M/j, which 
leads to particularly simple expressions, approximates Eq. p!7|) at the order of 1% for a 
reasonable range of masses [46] : 



1C\ 



5 17, 

- H h 

8 12 



(19) 



A particularly useful check of Eq. (jl7p is to evaluate it for q = t and expand it in the limit 
Mt > Mh-. 



M? cos a 



M? sin /3 



IC{ + 



My 



H . , MSUSY cos(a-^^, _ Ml ^.^^^ ^ ^ 
sin^ p Ml 



(20) 



1C\ 



/Co 



1C\ 



19 17, 1 

\ ks^ 

12 12 120Tt 



90rt ft 



+ 0{l/Tf) , 



1^ + 27^ + ^^^/^^^) 



(21) 



This result can also be directly obtained by asymptotic expansions of the corresponding 
Feynman diagrams along the lines of Refs. |31fl34j . Of course, we find complete agreement, 



Thanks to P. Slavich for confirmation. 
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therefore validating the result of Eq. (|17p. Let us stress that we do not use Eq. (|2ip in 
our numerical analysis, but rather the more general results obtained with the help of 
evalcsusy.f |32j . 

3 Total inclusive cross section 
3.1 Scenarios 

In this section we present numerical results for the hadronic cross section for Higgs pro- 
duction in gluon fusion in the MSSM. Through 0{ag), we include the effects from quark 
and squark loops as described above, i.e., the full quark mass dependence is kept for 
all leading order {a^q^) and real radiation contributions (Acjjj), as well as for the vir- 
tual quark loop terms (a^^^). The virtual top squark and top/stop/gluino contribution 
is taken into account in the limit {Mt, Mti, Mt2, Mg} ^ Mh, without further restrictions 
on the masses or SUSY parameters, with the help of evalcsusy.f. For the 2-loop vir- 
tual bottom squark and bottom/sbottom/gluino contribution, and only there, we apply 
the limit M^i = Mi,2 = Mg = Mf, ^ Mh,Mf), where in numerical evaluations we set 
Mfe = (Mfoi + Mb2 + Mg)/3. These approximations should be well-justified for the most 
popular SUSY benchmark scenarios (see, e.g. Refs. [30l[65]), in particular since they are 
only applied to the virtual effects where the partonic center-of-mass energy is fixed to 



In order to obtain numerical results for the cross section, one needs to insert numbers 
for the unknown SUSY parameters. In this paper, we will consider the following SUSY 
scenarios [301166]: 

gluophobic(±): 




MsusY = 350 GeV 
Ma = 300 GeV , 



/isusY =±300GeV 
M3 = 500GeV, Xt 



750 GeV 



(22) 



tan/3=10 



Mil ^ 150 GeV 
Mfci ^ 340 GeV 



Mt2 ^ 520 GeV 
Mb2 ~ 370 GeV 



Mg = 500 GeV 
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m. 



MsusY = ITeV. 
M2 = 200 GeV, 



tan/3=10 



Mti 



830 GeV 



/isusY = ±200GeV, 
M3 = 800 GeV , Xt 

Mt2 w 1170 GeV, 



2TeV 



(23) 



M, 



62 



1 TeV , M„ = 800 GeV 



no- mixing (lb): 



MsusY = 2TeV^ 
M2 = 200 GeV, 



//SUSY = ±200 GeV, 
M3 = 1600 GeV, Xt 



(24) 



tan/3=10 



i2 



b2 



2 TeV. 



M„ = 1600 GeV 



small Q!eff(ib): 



MsusY = 800 GeV , 
M2 = 500 GeV , 



/iSUSY =±2TeV, 



M, = 500 GeV 



Xt 



-1.1 TeV 



tan/3=10 



(25) 



Mti 
Mbi 



690 GeV , 
760 GeV , 



Mt2 : 
Mb2 



920 GeV , 
: 840 GeV 



M„ = 500 GeV 



The input parameters are the squark mass scale MsusY , the bihnear Higgs couphng fisvsY , 
the gaugino mass parameters AI2 and M3, and the off-diagonal term in the stop mixing 
matrix MfXt- Furthermore, it is always assumed that Ab = At ^ Xt+ hsusy / ^sm P , where 
Ab and At are trilinear couplings of the SUSY potential. For a more detailed description 
of these input parameters, let us refer to the documentation and references of the program 
FeynHiggs [67H70] which we use to determine the corresponding Higgs, sbottom and stop 
masses and mixing ang lesl The numbers for the squark masses quoted above are obtained 
for our default setting tan /3 = 10. For the quark masses, we use the input values 



Mt = 173.3 GeV , mb{mb) = 4.2 GeV Mb = 4.79 GeV , 

where Mt and Mb are on-shell masses, and nib is the MS running mass. 



(26) 



In the main part of this paper, we will restrict ourselves to the gluophobic(+) and/or the 
m™'"(+) scenario. The results for negative sign, including the other two scenarios, are 
deferred to Appendix ICl 

^Once the theoretical accuracy of the Higgs cross section increases further, the three-loop result for the 
Higgs mass has to be taken into account [43ll71lFf2j . 
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3.2 Numerical results 



The total inclusive Higgs production cross section within the MSSM through NLO QCD is 
shown for the gluophobic(+) and the ■mf^^^{+) scenario in Fig.|3](a) and (b), respectively. 
The figures include results for the Tevatron (pj5@1.96 TeV) as well as for the LHC at 
various energies, where tan /3 = 10. Here and in the following, the renormalization and 
factorization scales are set to ^ = = //r = Mh/2 [28], unless indicated otherwise. Also 
shown is the corresponding light Higgs mass M/^, see Fig.|3](c) and (d). The results for 
the other scenarios can be found in Appendix [Cl 

In the SM, where the cross section has been studied in great detail [9H22 t[23ti28l [73 l [7^. 
it was found that the NNLO QCD corrections are essential in order to reduce the scale 
uncertainty to an acceptable level. In addition, several calculations of beyond-NNLO 
effects lead to the conclusion that the fixed-order NNLO result provides a fairly precise 
prediction of the inclusive rate. 

Since we expect a similar behavior in the MSSM, we need to transfer the available infor- 
mation from the SM result to the MSSM case. We therefore define our best prediction of 
the total inclusive cross section as 



cr 



MSSM _ MSSM I ( h\2 



^MbbiVl I I „hVA /I , r \ SM,t SM,t /r)YN 



where c^lq^ is our result for the total inclusive cross section though 0{a\) within the 
MSSM as described in thepreceding sections. We consistently evaluate it with NLO parton 
density functions (PDFs) l2 The quantity c^j^^lo ^® top-quark induced SM cross section 
evaluated at (N)NLO QCD (i.e., with (N)NLO PDFs), and g^; = cos a/ sin /3. In addition, the 
electro- weak correction factor within the SM, (5ew [27j, is included by assuming complete 
factorization, as it was indicated to be a reasonable assumption [28]. 

In this paper we restrict most of our analysis to moderate and small tan /3. Therefore, we 
do not actually expect the resummation of the dominant tan/3 terms along the lines of 
Ref. |761I77] to be important. Nevertheless, we implement it, not least as a useful check: 
upon expansion of the resummed expression in terms of a^, we recover the coefficient K,\ 
(through 0{l/M^)) which gives the leading term in tan/3 at NLO. The numerical effect 
will be studied in more detail below. The renormalization and factorization scales are 
again set to //f = /ir = Mh/2. 

We provide the numerical results for 

^MSSM 

well as the individual contributions from 
Eq. ()27p for the gluophobic(-l-) and the m^''^^{+) scenario in Table [1] and [2j As expected, 
in the gluophobic scenario the cross section is typically much smaller than the SM value 



^In this paper, we use the central set of MSTW2008 [75] throughout. Detailed studies of the PDF 
dependence will be considered in a forthcoming paper. 
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gluophobic higgs(+) 




pp@ 14TeV 
pp @ 1 TeV 
pp @ 7 TeV 
pp@ 1.96 TeV 



1000 



pp@ 14 TeV 

m,"'='V) pp@10TeV 
" _pp@7TeV 
tanp = 10 pp@ 1.96 TeV 



100 



120 



140 160 
Ma [GeV] 



180 



200 



100 



120 



140 160 
Ma [GeV] 



180 



200 



(a) 



(b) 




100 120 140 160 180 200 100 120 140 160 180 200 

Mf, [GeV] [GeV] 

(c) (d) 

Figure 4: Inclusive total cross section for gluon fusion in the MSSM for the scenar- 
ios defined in Eqs. (j22p and (j23p at the Tevatron and the LHC for various energies, 
(a) gluophobic (+); (b) m™^^(-|-). Panels (c) and (d) show the corresponding 
light Higgs boson mass. 
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130 
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0.751 


1 
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1 
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1 
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155 


113.9 


0.902 


1 
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19.38 
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1 


.047 
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1 
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19.17 
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1 
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8.04 


175 
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1 
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180 
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1 


.047 
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0.960 


1 


.048 


16.23 


18.94 


4.81 


8.28 
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115.3 


0.964 


1 


.048 


16.19 


18.90 


4.87 


8.35 


195 


115.4 


0.968 


1 


.048 


16.17 


18.86 


4.92 


8.40 


200 


115.5 


0.971 


1 


.048 


16.14 


18.84 


4.97 


8.46 



Table 1: Data for the gluophobic(+) scenario with tan/3 = 10 at the LHC with 7TeV. 
(Masses are given in GeV, cross sections in pb.) 
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1 
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140 
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1 


.051 
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.052 


13.36 


15.54 


8.85 


11.31 


155 
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1 


.052 
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9.46 


12.08 
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0.913 


1 


.053 
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15.18 


9.69 


12.36 


170 
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0.928 


1 


.053 


13.00 


15.11 


9.88 


12.58 


175 


128.3 


0.939 


1 


.053 


12.95 


15.06 


10.03 


12.76 


180 


128.5 


0.948 


1 


.053 


12.91 


15.01 


10.16 


12.90 


185 


128.6 


0.955 


1 


.053 


12.88 


14.98 


10.27 


13.03 


190 


128.7 


0.961 


1 


.053 


12.86 


14.94 


10.36 


13.13 


195 


128.8 


0.965 


1 


.053 


12.83 


14.92 


10.44 


13.22 


200 


128.9 


0.969 


1 


.053 


12.81 


14.90 


10.51 


13.29 



Table 2: Data for the m^'"^{+) scenario with tan/3 = 10 at the LHC with 7TeV. (Masses 
are given in GeV, cross sections in pb.) 
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pp@10TeV pp@1.96TeV 




0.1 1 0.1 1 

(a) (b) 

Figure 5: Renormalization/factorization scale dependence {fi = fi{ = /^r) of the 
inclusive total cross section for gluon fusion in the MSSM. (a) LHC at lOTeV; (b) 
Tevatron. 



for the same Higgs mass, at least for Higgs masses above ~ 100 GeV. In particular for 
the phenomenologically most relevant region Mh ^ 114 GeV, the ratio of the NLO MSSM 
to SM result is only 25-30%. Inclusion of the top quark induced NNLO terms almost 
doubles the MSSM result, so that the MSSM to SM ratio increases to roughly 40-45%. For 
Mfi < 100 GeV, on the other hand, the MSSM cross section can become significantly larger 
than the SM result. 

This qualitative feature of a suppression of the cross section due to SUSY effects for 
Mfi > 114 GeV applies also for the m™'*^ scenario, but much less pronounced. And, also 
here, as one lowers the Higgs mass, the MSSM result surpasses the SM one. 

Fig. [5] shows the scale variation of the MSSM cross section in the m^^^{+) scenario at LO 
and NLO, as well as for our best prediction cj^^SSM (^labelled "NNLO"). Renormalization 
and factorization scale are identified in this plot (/i = ^Uf = /Ur), and varied by a factor 
10 around Mfi/2. The lower order results exhibit a similar behavior as known from the 
SM, and, as expected, the inclusion of the NNLO top quark terms leads to a considerable 
stabilization against the scale variation. The error due to scale variations estimated from 
this plot by considering the variation within M/j/4 < < Mh |28j is of the order of 15%. 

Fig. [6] shows the dependence of the cross section (Tevatron, 'm^^^{+)) on the choice of 
tan/3, where various contributions are displayed seperately. For example, the curve de- 
noted 6 + 6 is obtained by setting the top- and stop-Higgs couplings to zero, = g^^^ = 0, 
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pp@1.96TeV pp(3)1.96TeV 




5 10 15 20 25 30 35 40 6 8 10 12 14 

tan p tan p 

(a) (b) 

Figure 6: Dependence of (various contributions to) the cross section on tan /3 
for Tevatron conditions. Panel (b) is an zoom of the low- to intermediate-tan /3 
region of panel (a). 
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etc. The figure illustrates that already for tan/? ~ 10, the bottom effects outweigh the 
top effects (green/long vs. black/short dash-dotted). Furthermore, squark effects have a 
significant effect even at these large squark masses of 0{lTeV) if the dominant tan /3 terms 
are not resummed in the Higgs-bottom Yukawa coupling (orange/dash-double-dotted vs. 
purple/dashed). Once this resummation is taken into account, the squark effects become 
negligible, however (purple/dashed vs. red/solid). Note that this plot depends strongly 
on the choice of Ma, of course. Another well-known feature that is shown is the effect of 
the NLO corrections which are of the order of 100% (blue/dotted vs. red/solid). A very 
similar behavior is observed for LHC conditions, see Fig. [71 



4 Differential distributions 

In our perturbative partonic approach, the Higgs boson can only have non-zero transverse 
momentum when at least one parton is produced in association. Therefore, the purely 
virtual corrections do not contribute to the pT distribution da/dpT- In fact, the LO px 
distribution in the MSSM has been considered before |47H491[78|l79j . and we include it here 
only for the sake of completeness. In Fig.[8](a), the pT distribution at the LHC with 10 TeV 
is displayed for both the SM and the MSSM in the m™'^^(+) scenario Fig. [81(b) diplays 
the ratio of the two curves. We observe that the shape of the px distribution depends 
non-trivially on the model. The corresponding results for the Tevatron are shown in Fig.[9l 
The kink at pt ^ 150 GeV which is more pronounced in the SM, originates predominantly 

from the kinematical cut at \/I = \J Pj^ + + Pt (see also Ref. |53]). 
Concerning the rapidity distribution da/dy, where 

_ 1 E+p, 

with E and pz the energy and the longitudinal momentum of the Higgs boson in the lab 
frame, there clearly is a non-trivial distribution already at LO, given by 



da 
dy 



giVTey)giVTe-y) , r = ^, (29) 

LO ^ 



where s is the hadronic center-of-mass energy. Therefore, at NLO, also the virtual correc- 
tions need to be taken into account. We present here the first truly NLO results for this 
quantity in the MSSM. Fig.[ini(a) shows the rapidity distribution at the LHC with 10 TeV 
both for the SM and the MSSM in the m^'^^{+) scenario, while Fig.[TOl(b) shows again the 
ratio of the two curves. The model dependence of the shape is much smaller for the y- 



^Only the value of Mh influences the SM prediction when changing the SUSY parameters. 
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Figure 8: (a) Transverse momentum distribution of the Higgs boson at LO in the 
SM and the MSSM, using the m'^^{+) scenario at the LHC for lOTeV. (b) Ratio 
of the SM and the MSSM distribution. 



than for the pr-distribution. This is due to the fact that in the former case, a difference 
can arise only at NLO. 

The same conclusions hold for the Tevatron, for which the results are displayed in Fig.llll 



5 Fourth matter generation 



It is well-known that a SM-like 4* generation of quarks (denoted {t^, 64) in what follows) 
would increase the Higgs production cross section significantly, leading to a much larger 
exclusion region from the Tevatron search results [8]. On the other hand, the Higgs exclu- 
sion from electro- weak precision data would soften considerably with a 4*^ generation [80] . 
Effects of a 4*^ generation on the inclusive SM4 Higgs cross section have been evaluated 
through NNLO QCD [8HI82]. At LO, also H+iet and H+2 jet production has been studied 
in this model |82j . 

In this section, we evaluate the effect of an MSSM-like 4**^ generation on the Higgs produc- 
tion cross section through NLO. Of course, also the SUSY constraints on the Higgs mass 
change significantly with the presence of a 4**^ generation |83H85j . For this first study, we 
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Figure 9: (a) Transverse momentum distribution of the Higgs boson at LO in the 
SM and the MSSM, using the m™^^(+) scenario at the Tevatron. (b) Ratio of the 
SM and the MSSM distribution. 




Figure 10: (a) Rapidity distribution of the Higgs boson at NLO in the SM and 
the MSSM, using the m^^^(+) scenario at the LHC for 10 TeV. (b) Ratio of the 
MSSM and the SM distribution. 
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pp @ 1 .96 TeV pp @ 1 .96 TeV 




(a) (b) 

Figure 11: (a) Rapidity distribution of the Higgs boson at NLO in the SM and 
the MSSM, using the m™^''(+) scenario at the Tevatron. (b) Ratio of the MSSM 
and the SM distribution. 



will use the approximation 

Mi = iMt'^^f+ E M| = M|usY + M|, (30) 

q=t4,b4 ^ 

where v = 246 GeV, MsusY is given by the scenarios as defined in Sect. 13. ll and M4 is the 
mass of the fourth generation quarks which are taken to be degenerate. The 3-generation 
result M^^^^ is again taken from FeynHiggs [67fl70j . Furthermore, for the two lighter 
4th generation squark masses we also assume Mgi = M4 (g G {64, t4}), while the heavier 
ones are set to Mq2 = (2MgugY + M^)^/^. The form of the coupling constants for the 4**^ 
generation (s)quarks is taken to be identical to the first three generations, i.e., they are 
given by Eqs. (f33]l and (f37|) of Appendix lAj with the obvious replacements of the masses 
and mixing angles. For recent 4*^ generation search limits, see Refs. j86y87j. 

Fig.[T2](a) displays the Tevatron and LHC results for the SM4 and the MSSM4 using the 
SUSY parameters as in the ?n™'^^(+) scenario. In order to avoid non-perturbativity of 
the /16464 coupling, we set tan/3 = 1. On top of the well-known factor of ~ 9 from the 
SM— )-SM4 transition, we find a 30% decrease when switching to the MSSM4. This value 
depends quite sensitively on the actual values of the 4*^ generation squark masses and 
suggests further, more detailed investigations. 
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Figure 12: Higgs production cross section as a function of the light Higgs boson 
mass (a) at the Tevatron and (b) at the LHC, with a 4**^ SM- and MSSM-hke 
generation of (s)fermions. 



6 Conclusions 

The various contributions of the NLO QCD corrections to the cross section for Higgs 
production in gluon fusion have been combined in order to evaluate it consistently within 
the MSSM. Both quark and squark effects have been taken into account in the top and 
the bottom sector, including the mixed quark/squark/gluino and all interference effects. 
The numerical results have been presented for a set of selected MSSM scenarios, but our 
implementation allows for any other reasonable set of parameters. 

For the total inclusive cross section, we constructed a "best approximation" which includes 
the known NNLO QCD top-quark induced corrections. We also provided results for Higgs 
distributions in px through LO and in y through NLO. Finally, the effect of a 4**^ MSSM 
generation on the Higgs production rate has been investigated. 

In future works, we will study the effect of changing the renormalization scheme, in par- 
ticular for the bottom/sbottom sector, include the recently evaluated NNLO effects for the 
top/stop sector, and consider also the production of the heavy Higgs (see also Ref. |38j). 

We hope that these results, plus similar ones provided on-line [88], will be useful input 
for experimental analyses. We will gladly provide numbers for other scenarios to the 
interested reader, and are planning to release the code for public use in the near future. 
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A Feynman rules and coupling constants 



The Feynman rules for the SUSY-QCD vertices can be found in Refs. |32p61j. Here, we 
give in addition the ones for the bottom- and sbottom-Higgs vertices: 




The bottom-Higgs coupling constant reads 



h sma , , 



and the sbottom-Higgs couplings are 

h _ h,EW , h,ii . h,a 



9b,ij = 9b,ij + 9b fj + 9b,ij > (32) 
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with 



h,EW 

9b,n 

/i,EW 
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/i,EW 
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h,ii 
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%,22 
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/i,EW 1 / 
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cff) Sin 206. 
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"56 22 = ^ ---on — sin 26'6 , 
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56,21 



M6 cos2 /3 
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M6 



COS 2^6 ) 



sma 



cos j3 



sma 



2 + 



61 



2M2 



^ sin^ 26^ 



cos ;S 

h,a 
56,21 



^6^1 



2M2 



M.2 „ 

^ sin^ 2^5 



sinaM,\-M2, ^.^^^^^^^^^^^ 



cos /3 



2M2 



(33) 



where 



9 M'^ 

c!,r = ^ (|/3| - IQI sin^ 0w) sm{a + /3) , 
9 A/f2 

c^^ = -rjflQl sin' 9w Ma + (34) 



sin0^ = ^/l-^, /3 = -^, Q = -^ 



and 



tan/3 = ^; = ^^ = ^^= = y^^27^^ 246GeV, (35) 

with fi, the vacuum expectation values of the two Higgs doublets. The angle ^6 rotates 
the super-partners of the left- and right-handed bottom quarks into their mass eigenstates, 
while a does the same for the neutral components of the Higgs doublets. In Eq. (I33p we 
have already expressed the trilinear couplings of the soft SUSY breaking terms through 
independent parameters: 

Ab = sin 29, + /xsusY tan /3 . (36) 
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The expressions in Eq. (j33p and (j34p are completely analogous to the ones for the scalar top 
sector, given in Ref. [32]. They can be transformed into each other by the replacements 



t^b, /^SUSY ^ -AtsusY , 

Q o Q - 1 , /3 ^ /3 - 1 . 



vr 

a o a + — , 



vr 



(37) 



B Expansion of quark terms 



Analytic expressions for the virtual quark loop 
at NLO have been presented ir|l Refs [MIIMIE 
have 
limit, 



Lc viiULiai v^aark loop contribution to the gluon fusion process 
ah iNuvv have been presented irH Refs [35 1 136 1 144]. It may be convenient for the reader to 
have the relevant limits of these expressions. We quote them here for the large quark Tnn«« 



mass 



Uq _ 
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— = — h 

,h A ^ 



11 1237 



1 17863 2 
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122835309478756875 " 316945428161236875 " ^ " ^ 



and for the small quark mass limit: 



9^ 



= Tl 
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+ L 



47 - 31 Cs - 18 C4 + (-12 - 4C2 + 11 Ca) + (-^ - \ C2) 



48 
11 



1 - 10 C2 + 128 Ca + 72 C4 + (-2 - 44 Cs) 



y + 2C2) 
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r. 



q) ) 



+ 0, 

where = C(^) is Riemann's zeta function evaluated at n, with 



it" 

(2 = — = 1.64493 . . . , Cs = 1.20206 
6 



C4 = ^ = 1.08232 . 



(38) 



(39) 



(40) 



*For the sake of clarity, let us remark that in the notation of Ref. 



it isFo^Bf = 2aW/g^ 
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C Total cross section for other SUSY scenarios 



This section collects the results for the total inclusive cross section in the SUSY scenarios 
of Eqs. (j22p - (j25p which have not been discussed in the main text. We restrict ourselves 
to graphical representations; the corresponding data tables can be found at the URL |88j . 
Figs. ll3H15l are simply the analogues of Fig. U] for these scenarios. 
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found in Fig. HI 
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Figure 14: Inclusive total cross section for gluon fusion in the MSSM. (a) no- 
mixing(— ); (b) no-mixing(+). Panels (c) and (d) show the corresponding light 
Higgs boson mass. 
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Figure 15: Inclusive total cross section for gluon fusion in the MSSM. (a) small- 
Oicsi—)', (b) small-Qcfr(+). Panels (c) and (d) show the corresponding hght 
Higgs boson mass. 
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